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Elements of length, area 

Exertia, («i/? — 1) 

Force per unit volume of a phase 
Body force per unit mass 
Unit dyadic 
Volumetric flux density 
Kinetic energy density 
Geurst’s added mass coefficient 
Momentum density 
Added mass force 
Mean pressure in a phase 
Macroscopic pressure 

Combined momentum flux and stress tensor 

Density ratio P 1 /P 2 

Time 

Microscopic velocity 
Average velocity 
Relative velocity Vi — V 2 
Volumetric fraction of a phase 
“Resistivity,” (1) 

Macroscopic potential 
Specific momentum, (24), (25) 

Density 

Macroscopic potential 
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Subscripts 


With particles at rest 
Phase 1 
Phase 2 
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INTRODUCTION 


Potential flow theory for a single fluid has been established for many years. Although 
its limitations for describing real motions are well known, it does provide a self-consistent 
structure for analysis and often provides a reasonably accurate description of at least part 
of the flow field of fluids with low viscosity. 

As an example of two-phase flow, one can imagine a suspension of particles in a fluid 
that obeys all the requirements of potential theory at the microscopic level. If the only 
forces acting on the particles are “conservative,” it would appear that their motion might 
also reasonably be expected to be describable in terms of a suitable potential. Averaging 
of these potentials would lead to macroscopic potentials, true properties of the mixture, 
that should be related in some way to the average motion of the phases. Indeed, previous 
attempts to determine the inertial coupling terms in the two-fluid model have implicitly 
assumed potential flow at the microscopic level. 

A complete two-phase theory of this type will be as idealized as was the classical 
theory of single-phase potential flow. However, it should be useful in the same sorts of 
ways, both as an approximation in many situations and as a standard that must at some 
level be consistent with other approaches, such as those that attempt to define closure 
relations for a set of averaged basic equations. 

This paper describes some features of two recent approaches along these lines. The first 
is based on a set of progressive examples that can be analyzed using common techniques, 
such as conservation laws, and taken together appear to lead in the direction of a general 
theory, the tactic used in [1], The second is based on variational methods , a classical 
approach to conservative mechanical systems that has a respectable history of application 
to single phase flows. The latter approach, exemplified by several recent papers by Geurst 
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[2-5], appears generally to be consistent with the former, at least in those cases for which 
it has been possible to obtain comparable results. 

PROGRESSIVE EXAMPLES 

The theory developed in [1] starts with a situation where the particles are at rest. 
Fluid flows past these particles as though a porous medium. At the microscopic level 
the equipotentials are not “smooth” but they do not differ much from the more gross 
macroscopic equipotential. Differences between the two levels of equipotential are smeared 
out over lengths comparable with the particle size. I do not have a rigorous proof, but I 
believe it is valid to treat these two equipotential surfaces as essentially identical for most 
purposes. This is not true of the other properties, such as pressure and velocity, that vary 
more at the microscopic than the macroscopic scale and must be averaged carefully. 

Just as in electrical conduction past a matrix of non-conducting particles, the macro- 
scopic fluid flux will be proportional to the macroscopic potential gradient, the “resistivity” 
being represented by a factor 0 that depends only on the particle geometry and the void 
fraction, as long as the arrangement is isotropic. We therefore have 

jo = (1) 

Since 0 is unity for unimpeded flow, it will be greater than one when particles are 
present. A crude description of the situation could be to say that some of the fluid is “held 
up” or “entrained” by the particles so that only a fraction of the space is available for 

direct flow. 

The average velocity of the fluid is the relative velocity , 

w = ^ = (2) 

Oi 0. 10 
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The kinetic energy of the fluid in unit volume confined between two equipotentials is, 
from a standard theorem of potential flow. 



/ 


<^u • ds = • jo = ]^p\ct\Pw 2 


The kinetic energy per unit volume of phase 1 is 


(3) 


ki 



- i piw 2 (a 1 ft) 


(4) 


If we denote the fluid velocity at the microscopic level by m , (2) and (4) are equivalent 


to 


W =< Uj > 


(5) 


w 2 (*ip =< Uj > (6) 

Invoking the Schwarz Inequality, it is clear from (5) and (6) that cn/3 is greater than 1, a 
proof pointed out to me by my student Chao Luo. 

By considering the changes in the fluid kinetic energy resulting from a uniform volume 
change for every particle it can be shown [1] that the difference is mean pressure between 
the phases is 


P2 — Pi — ~P\w 2 a\ ~ 1 ) 


(7) 


The additional “1” in (7) appears to be gratuitous. It is introduced at this stage 
because the “exertia,” defined [1] as 
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£ = (<*!/?-!) 


( 8 ) 


turns out to reappear in numerous contexts and to be one version of an added mass 
coefficient. 

If we now superimpose a uniform velocity V 2 on the above motion, this is equivalent 
to superimposing a corresponding potential gradient and we obtain 


Vi + Ew = — 


(9) 


where the relative velocity is 


w = Vi — v 2 


( 10 ) 


The net momentum and kinetic energy densities are then 


m — p\ 0 c\\\ 4* P 2 & 2 V 2 


(ii) 


k = 


^PlOilvl + ^ piOtiEw 2 


( 12 ) 


From (11) and (12) we may deduce that the effective equations of motion of a uniform 
suspension accelerating under the influence of a uniform macroscopic pressure gradient 
and body force fields are 


25 


V] + Ew — 


VP 


ai xw 


(13) 


Pi 


+ gi + 


Pi 


Pl<*l „ • VP , , a 2><W 

v 2 — —E w = — 1- g 2 + 


p2&2 


P 2 


P2 


ai and a 2 are arbitrary vectors subject to the constraint 


(14) 


(c*iai + o 2 a 2 )xw = 0 (15) 

It would appear that the assembly must be incompressible if it is to be truly “uniform” 
in a pressure gradient. If gi is a conservative force field, comparison between (9) and (13) 
would seem to indicate that Vxai x w is zero and therefore Vxa 2 x w is as well from (15). 
Eq. (14) then suggests that if g 2 is conservative, the combination on the left-hand side is 
the gradient of another “potential” which we could define in a form somewhat resembling 
(9): 


v 2 - : — V»7 (16) 

P 2 £*2 

A different argument is used to derive the one-dimensional equivalent of (16) in [1]. In each 
case, the interpretation of $ and r) could be, in the classical view, in terms of impulsive 
pressures and body forces necessary to set up the motion. This derivation, made for a 
uniform suspension, requires further argument, or a leap of faith, if it is to be applied 
more generally to a dispersion in which w and oi, and hence /? or E, vary with position. 

Eqs. (13) and (14) show that the exertia plays the classical role of a coefficient of 
apparent mass , being proportional to various alternative definitions [1] in the literature. 
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Some other results derived from mechanistic arguments [1] axe: 

• The combined momentum flux and stress tensor for the suspension: 


p = «i(pi V 1 V 1 + pi I) + a 2 (p2V 2 V2 + P2 1) + aipiEw-w (17) 

• Bernoulli’s equation for fluid flowing steadily past a stationary particle matrix: 

Pi = Poi - ^P\w 2 {\ + E) (18) 

• The force per unit volume on a particle in a stationary lattice: 

f 2 = Vp 2 = V ( 19 ) 

When these results are used to check several hypothesized forms of the equations of 
motion in the two- fluid model of two-phase flow, there axe found to be discrepancies [7], 
except when Geurst’s equations are used. 

GEURST’S EQUATIONS 

In a series of papers [2-6] Geurst has used variational methods to derive a number of 
results, starting from the hypothesis that the kinetic energy per unit total volume is 

k = ^piaivl + ^p 2 a 2 vl + ^pimw 2 ( 20 ) 

which is the same as (12) with the alternative definition 


m = ai(c*i /? — 1) = a\E 


( 21 ) 
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Initially Geurst [2,3] developed a one- dimensional theory, introducing Lagrange mul- 
tipliers that are the same as the potentials $ and rj in (9) and (16). He went on to 
derive equations of motion, momentum, fluxes, pressures, etc. that are compatible with 
the results in Section 2, where there is a clear equivalence. 

In three-dimensions the derivation is complicated by the introduction of Clebsch po- 
tentials and Lin constraints. Details are not provided in earlier publications [2-4] but they 
do appear in a recent one [6] in which the equivalents of (9) and (16) are expressed as 


TTi = V<f> 1 + tfiVxi (22) 

*2 = V<f>2 + 02^X2 H V<^ n (23) 

P2 

where 7Ti and 7T2 are generalized specific momenta defined as 

ui 

Tfi = Vi (v 2 -vi) (24) 

Otl 

P\m v , . 

7T2=V 2 + - (V 2 -V X ) (25) 

p2&2 

which are identical with the left-hand sides of (9) and (16), in view of (21). 

The <^’s and x’ s are fi ye “potentials” or Lagrange multipliers and n is the bubble 
density. The detailed interpretation of these terms is not important in the present context. 

Geurst’s equations of motion appear in many equivalent forms. The versions selected 
in [8] may be written in the present notation as 
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— (piaiwi) + V • (piaiYiKi) = -«iVp 2 - V 


' V ...2 


~pi(m + aim )w 


+«ifi + Pim (\2 - vi) • (Vv 2 ) J 


( 26 ) 


— (p 2 «27T2) + v • (p 2 a 2 v 2 ir 2 ) = ~a 2 Vp 2 + a 2 f 2 - p\m(\ 2 - Vj) • (Vv 2 ) T (27) 
at 

which have the form of conservation laws for and 7T2. m* denotes dm/da.^. 

The two continuity equations are 

$ 

— (/>iori) + V • (/9iaiv x ) = 0 (28) 

d 

— (p 2 a 2 ) + X/-(p 2 a 2 v 2 ) = 0 (29) 

If (7) is accepted and (21) is used, we find that the term involving w 2 on the right-hand 
side of (26) may be rewritten as 

\p\{jn + ctira')u> 2 = ai(pi - p 2 ) (30) 

Using (28) through (30) it is straightforward to arrange (26) and (27) into the forms 


d_ 

dt 


(7Ti) + V ( Vi 


71-1 2 


\~{^2 -VO 2 ) Vj X V X7Ti 


Vpi + £l 

p\ pi 


(31) 


d_ 

dt 


(tt 2 ) + V 


^V 2 • 7r 2 



v 2 x V x ir 2 


— ^P2 h_ 

P2 p2 


(32) 
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Eqs. (31) and (32) are versions of the Bernoulli Equations derived in [8] without 
identifying 7Ti and rc-i with gradients of potentials. 

Now, if either p\ is constant or Vpi is parallel to Vpi, and V x {i\j p\) = 0, we may 
take the curl of (31) and obtain 

— (VxiTi) — Vxvx x(VX7Ti) = 0 (33) 

ot 

which is a conservation law for the vorticity of 7Ti . If we consider the rate of change of the 
net vorticity of threading a loop moving with the velocity Vi , we obtain 

s/’ r ‘ d£ =l/ (Vx ’ ri)ds 

= j ^(VxTT]) • ds 4- J (Vxjti)- (vi xdf) (34) 

The two terms in (34) represent the sum of in-place changes in the flux through the loop 
and contributions picked up by the motion of boundaries. Changing the order of the scalar 
triple product in the final term and invoking Stokes’ Theorem, we obtain a surface integral 
of (33) which is identically zero, in other words, 

j t J 7n • di = 0 (35) 

This result was obtained by Geurst [3]. A similar conclusion follows for (Vx-n^). 

Now, if both phases come from a region in which there is no curl to 7Ti and 7T2 (for 
example, a stagnation region), conservative body forces act and density gradients (if any) 
are parallel to pressure gradients, then throughout the flow, in view of (35), it should be 
true that 
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Vx7Ti = 0; 


V X7T 2 = 0 


(36) 


which implies that both 7Ti and 7T2 are gradients of suitable potentials and only the first 
term is needed on the right-hand sides of (22) and (23) which reduce to the form of (9) and 
(16). This development apparently makes more explicit the conditions for the existence of 
two-phase potential flow , requirements which parallel those for the classical single-phase 
case. 

DISCUSSION 

The previous sections of this paper have outlined two approaches to potential two- 
phase flow. Each has a justifiable theoretical base and is self-consistent. Moreover, both 
approaches appear to give the “right” prediction for several well-defined situations [7] while 
some other formulations fail these tests. 

In order for these ideas to blossom further there need to be: 

a) Further developments, from the same basic set of assumptions, that encompass 
more generality. 

b) More rigorous derivations that clearly explain the order of approximation involved 
in treating the flow of discrete entities as a continuum. 

c) Reconciliation with alternative approaches, particularly those involving averag- 
ing. 

d) More solutions to specific problems that can be thoroughly investigated for con- 
sistency. 

e) An understanding of outstanding incompatibilities between these approaches and 
various other theories, with a clear explanation of what has “gone wrong.” 
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I hope that some of these items can be discussed at this meeting. At this time I only 
wish to briefly address (d) and (e). 

Apart from the “tests” described in [1] and [7] and perturbation techniques leading 
to the description of wave propagation and stability in bubbly flow [2,3,5], I know of 
only two complete solutions to Geurst’s equations. The first involves steady flow of two 
incompressible fluids from a common stagnation region. The simple conclusion that is 
reached is that the void fraction is constant and the velocities of both phases are the 
gradients of potentials that are proportional to each other and may be borrowed from an 
equivalent single-phase flow. The velocity ratio, or slip ratio is 


where 



(37) 


R = pi I pi 


(38) 


In the above derivation use was made of Maxwell’s [9] approximate expression for the 
exertia , 


Es3 Y ( 39 ) 

The second solution again uses (39) and leads to a relationship between the velocities 
in unsteady incompressible flow [8] 
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This has the same form as a relationship which would be obtained by eliminating the 
pressures from equations of motion that ignored inertial coupling, except that the effective 
densities are changed. There is opportunity to test (40) by comparison with the one- 
dimensional transient response of fluidized beds. 

Regarding item (e) I will repeat, in a slightly different form, Geurst's equations of 
motion [4] as presented in [8] 


— (piaiVi) + V • (pioivivi + pimww) — 
dt 

-aiVpi - (pi - p 2 )Vc*i + Mf + aifi 


(41) 


— (p 2 « 2 v 2 ) + V • (p 2 a 2 v 2 v 2 ) = -« 2 Vp2 - M ? + a 2 f 2 
dt 

with the added mass term expressed as 


(42) 


M( = p\m 


(| +V 2 . V )v 2 -(|+ Vi v ) v 1 + 


wxVxvi' 




+p\mV— + w 


— (pim) + V • p\m\2 
dt 


] 


(43) 


This set of equations contains several more terms than one would find in most similar 
expressions in the literature. By dint of these extra terms, several “tests" are passed that 
other formulations fail (7). Specifically, these tests involve (18) and (19), which reduces 
for a dilute dispersion to Taylor’s expression |10] for the force on a stationary object in 
an accelerating flow. It would be desirable to devise other “tests” that might help to 
discriminate further between true and false expressions. 
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